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Electronic energies and wave functions of ground and excited states of LiO are calculated in the 
Hartree-Fock approximation The vibrational oscillator strengths of the X 2 I1 and A 2 S + valence states 
are calculated for v= to 1 and 1 to 2 transitions. Electronic oscillator strengths are also presented for 
transitions to the quantum number two Li Rydbergs. Reasons are presented to support the use of the 
Hartree-Fock approximation. 
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1 . Introduction 

Refractory metal oxides have proved to be difficult 
experimental subjects. For even one of the simplest, 
LiO, little is known about the system beyond the ground 
state properties of the X 2 I1 state. In particular there 
is no experimental information on the infrared intensi- 
ties and the ultraviolet spectrum. Interest in the radia- 
tive properties of LiO has prompted a number of 
theoretical calculations. Accurate Hartree-Fock (H-F) 
and configuration interaction (CI) results [l] 1 have 
been reported for the X 2 H and the low-lying 2 £+ excited 
states. These results show that the properties of this 
molecule are accessible to theoretical calculation. This 
paper will apply the Hartree-Fock technique to some 
remaining problems of the properties of LiO. 

The Hartree-Fock self-consistent-field variational 
calculations are based upon an anti-symmetrized 
product of one-electron orbitals that are determined by 
the averaged field of all the other electrons. The lack 
of correlation in this model leads to large absolute 
errors in the energy. But the applicability of the Hartree- 
Fock technique must not be rejected out of hand with- 
out an examination of the problem at question since 
correlated wave functions are obtained at high cost. 
There are four calculations considered in this paper: 

(1) the energy curve of the ground state and the verti- 
cal excitation energies of the low-lying excited states, 

(2) the infrared intensities of the X 2 n and 2 X + states, 

(3) the oscillator strengths for excitation from the ground 
state to all the low-lying excited states, and (4) the A 
doubling matrix elements. 



The limitations of the Hartree-Fock method can lead 
to significant and unknown errors in the calculation 
of excitation energies when the correlation energies 
are very different in the ground and excited states. 
However, it is clear from the monovalent character 
of the Li and the ionic behavior of the molecule in the 
neighborhood of the equilibrium internuclear separa- 
tion, that the X 2 n and A 2 2 + states are best described 
as different orientations of the 0~( 2 P) orbital angular 
momentum relative to the molecular axis. Hartree- 
Fock and correlated results should agree for relative 
energies and this is found in the recent accurate CI 
results [lb]. The present study supplements this work 
done on the X 2 n and A 2 £+ states and repeats sufficient 
work to show substantial agreement in the areas that 
overlap. 



2. Self-Consistent-Field Solutions 



The present calculations were performed with a 
Gaussian-type function (GTF) constructed by aug- 
menting the 9/5 basis of Huzinaga [3] with three 2p Li 
functions that are suitable for representing the 2 P 
excited state of Li. The configurations can be repre- 
sented: 



lo- 2 2o- 2 3o- 2 4o- 2 l7r 3 X 2 n 



lo- 2 2o- 2 3cr 2 4o- Itt 4 A 2 2+ 



* This research was supported in part by the Advanced Research Projects Agency of the 
Department of Defense under the Strategic Technology Office. 

1 Figures in brackets indicate the literature references at the end of this paper. 



The X 2 n and A 2 2 + were calculated at five distances 
for which the energies and dipole moments are given 
in table 1. There is a crossing of these states near 
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Table 1. LiO Energy as a function of internuclear separation* 





x 2 n 




A*£+ 




/?(a.u.) 


-E 


M 


-E 


M 


2.7 


82.2493 


2.49 


82.2586 


2.19 


2.9 


82.2706 


2.63 


82.2699 


2.37 


3.184 


82.2817 


2.84 


82.2716 


2.64 


3.4 


82.2815 


2.99 


82.2668 


2.85 


3.6 


82.2774 


3.14 


82.2598 


3.04 



a Energy and dipole moment given in atomic units: E, 1 
a.u. = 27.2097 eV; /i, 1 a.u. = 2.54158. KH 8 esu cm. 

2.9 a.u. which is in accord with the united atom char- 
acter of these states. At large internuclear separation 
the 2 2+ state correlates with 0( 3 P) + Li( 2 P) while 
the 2 I1 state correlates with ground state atoms. Since 
the SCF configurations do not go correctly asymp- 
totically there is little point in determining dissociation 
energies. 

The Rydberg orbitals were obtained by direct 
variation at /? = 3.184 a.u. of the Rydberg states to the 
3 S~ and J X + ions. These states were chosen as repre- 
sentative by consideration of the ion atom interaction 
at long distance. Adiabatic correlation rules predict 
that the interaction of Li+( X S) and 0( 3 P) yields 3 S~ 
and 3 I1 states for which the 3 £~ is attractive and 3 1I 
repulsive with respect to the ion-quadrupole inter- 
action. Similarly, the interaction of Li + ( J S) and 
O0D) yields *A, i n, and *£+ for which only the *A 
state is attractive with respect to the ion-quadrupole 
interaction. The most attractive ion states were con- 
sidered with 2s and 2p Rydberg Li electrons as given 
in the list above. Analysis of the wave functions shows 
very little difference between Rydberg orbitals obtained 
between the Rydbergs for the 3 2 _ and x £+ states. 
The different orientation of the quadrupole of the 
oxygen atom will affect the energy position but the 
transition probabilities for a number of Rydberg 
transactions can be estimated from the orbitals calcu- 
lated from a few representative states. 

3. Discussion 

A. X 2 IIcind A 2 2+ Energies 

Polarization functions were not included in the 
present calculation. Comparison with the accurate 
Hartree-Fock results of Yoshimine [lb] show that the 
error increases to smaller internuclear separation as 
expected. Polarization is more important as the 
atoms approach each other. However, the error induced 
in the spectroscopic properties is not substantial. 
This calculation obtained o) e = 804.4 cm" 1 , R e = 3.275 
a.u. , and o) e x e = 6.1 cm -1 for X 2 II and o) e = 877.5 cm -1 , 
R e = 3.076 a.u., and co e x e = 10.5 cm" 1 for A 2 £+. The 
largest error is about 5 percent in the (o e of X 2 II. 
Note that the accurate CI spectroscopic constants 
are in very close agreement with the H-F results [lb]. 
Fortuitously, the present results are in better agreement 



with the reported experimental results [2] for the 
X 2 I1 state. The (o e for 2 2 + is within 7 cm -1 of the 
H-F value and the improved shape is related both 
to the fact that the asymptote of the 2 S + H-F state 
is not ionic which diminishes the asymptotic error 
and to the cr character of the open shell orbital. The 
GTF basis is more flexible for the cr orbitals and the 
polarization error is accordingly less. 

B. Infrared Intensities 

The dipole moment values are within 10 percent 
at worst in absolute agreement with the accurate CI 
values [1]. More important all the curves are very 
ionic and linear. The slopes of the present results 
agree with the CI results to within 2 to 3 percent. 

Vibrational oscillator strengths were calculated 
from the equation 



/=-A£|<¥, 



= 3 A£/2 ' 



|t(Jt)|«V>| 



The vibrational energies and /numbers for the 0— ► 1 
and 1 — » 2 transitions are given in table 2. The infrared 

Table 2. LiO vibrational oscillator strength 



State 



Vibrational 
transition 



kE{cm- 1 ) 



Band 

strength 

<P) 



MO 5 



x 2 n 


0-1 

1-2 


792 
780 


0.0082 
.0165 


2.0 

3.9 


A2£ + 


0-1 
1-2 


856 
836 


.0128 
.0260 


3.3 

7.4 



transitions are strong for both states as expected for 
an ionic system. These values are in moderate agree- 
ment [4] with values obtained from the CI calculation. 
Since the dipole moment function is linear the 0— ► 2 
transition is very weak. 

The vibration-rotational levels of the X 2 I1 state are 
perturbed by the A 2 2+. The A-doubling matrix element 
is readily evaluated since the electronic matrix 
element <^(7r) \ L + |¥(2+) > does not vary with 
internuclear separation and is equal in magnitude to 
the pure precessing value. The matrix element 
<¥„( 2 II) |F|'«V(«2+)>, where 



V = 



1 



ViixR 2 



is given in table 3. The sensitivity of these values with 
vibrational quantum number would caution that they 
be considered quite approximate considering the level 
of accuracy of the energy curves. 
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Table 3. LiO A-doubling matrix element '' 



v( 2 U) t/( 2 2 + ) 





1 


2 


3 


4 





1.25 


1.03 


0.69 


0.43 


0.25 


1 


-0.96 


0.47 


.90 


.83 


.62 


2 


0.40 


-0.99 


-0.008 


.65 


.81 


3 


-0.09 


0.57 


-.89 


-0.29 


.39 


4 


0.0007 


-0.14 


0.67 


-.77 


-0.44 



■ Matrix element < iM 2 n) \h 2 l V2)u,R 2 |i/v ( 2 2 + ) > in cm- 1 . 

C. Electronic Oscillator Strengths 

The low-lying ion states of LiO+ can be represented: 



1o- 2 2o- 2 3o- 2 4o- 2 1tt 2 



*A, ^ 



lo*2o*3o*4crl'n 3 3 II, *!! 



Rydberg Li 2s and 2p orbitals were generated in the 
field of the 3 £~ and 1 S + SCF approximations to the 
ions. Since these states are most widely separated 
in energy and the Rydberg orbitals were determined 
to be substantially the same for both ion cores, the 
Rydberg 5o-(Li 2s), 6o-(Li 2p), and 27r(Li 2p) orbitals 
were used to determine oscillator strengths for the 
3 2~, J A, 3 II, and *n Rydbergs. The X S + state will 
require substantial mixing of the 4cr 2 l7r 2 and Itt 4 
configurations and, in any case, is the highest energy 
of these ion states. 

Asymptotic correlations are given in table 4 for 



Table 4. Adiabatic correlation of molecular states of LiO 



Molecular states a 


Atomic states b 


Asymptotic energy 

(eV) 


*< 4 2+,2-(2), 


Li(^P„)+0(»P 9 ) 


3.83 


n(2),A 






2 ' 4 2-, n 


L1(»S,)+0('P,) 


3.37 


2 n(3),s-, x + (2), <f> 


Li(2P„)+0('D g ) 


2.57 


A(2),<& 






2 A, s+, n 


U(»S,) +0(»D,) 


1.97 


2 ' 4 S-(2),n(2), 


Li(2P„)+0(3P 9 ) 


1.85 


A, 2+ 






2 ' 4 2-, n 


Li(«S,)+0('P,) 






a2 > 4 S-, n signifies *2-, 4 2", 2 I1, 4 ll. 

b Atomic states and energies, from C. E. Moore, Atomic Energy 
Levels, Nat. Stand. Ref. Data Ser., Nat. Bur. Stand. (U.S.), 35, Vol. 1, 
359 pages (Dec. 1971). 



several of the lower asymptotes that can correlate 
with the valence and quantum number 2 Rydbergs. 
It is apparent that the asymptotic behavior of the 
Rydberg states is complicated since a sizable fraction 
must become valence. However, at the R e of the ground 
state the single Rydberg configuration should be a 
good approximation to the true wave function. 

From the asymptotic behavior and the calculated 
energies at R = 3.184 a.u., we can deduce that these 
Rydberg states will probably be photodissociated. 
Both the 4 S~ and 2 £~ that arise from the ground state 
atoms will correlate with the l7r 2 5cr configuration and 
the 4 S~ state is at least 0.6 eV above the asymptote. 
The 2 A state correlates to Li( 2 P) + O ( 3 P) and is calcu- 
lated to be about 1.0 eV above this asymptote. The 
ion states go formally correctly to the asymptote, and 
the configuration mixing required to insure that the 
excited valence asymptotes are correctly obtained 
should only be important at large internuclear separa- 
tion. The excitation energy is, therefore, meaningful 
within the error due to the lack of polarization func- 
tions. The error in the X 2 II state at R e is about 0.027 
a.u. and the 2 A state would still be repulsive at this 
distance. The 2 2~ and 2 II Rydbergs that can go asymp- 
totically to Li( 2 P)-hO( 3 P) should be above the 4 2~ 
in energy and, even with the polarization error, would 
be photodissociated. Without a configuration mixing 
calculation it is impossible to be certain but the present 
analysis would indicate that the vertical excitation of 
LiO from either X 2 I1 or A 2 2+ would lead to dissociation 
throughout the Rydberg region. 

Transition probabilities are calculated from both 
the X 2 n and A 2 2+ states by forming the appropriate 
symmetry-adapted linear combination of determinants. 
In order to avoid estimating the excitation energy, 
the oscillator strengths were obtained with an ex- 
pression using the product of the dipole length and 
velocity matrix elements [5], 



/=3 |<r> 2 <p> 2 |i/ 2 (2- 



■So.\ i 8oAy)/(2- 



So A-) 



after summing over all initial and final state degen- 
eracies, averaging over the degeneracies of the initial 
state, including the A-doublet substates [6], and using 
the fact that the sum of Franck-Condon factors to 
final vibrational states is unity. 

There are no strong transitions apparent in table 5. 
The strong 2s— »2p resonance transition of Li is not 
involved here since the ground state is so ionic that 
all low-lying transitions originate on the 0~ atom and 
the overlap through the dipole operator is relatively 
small. 

The electronic oscillator strength for the transition 
from the ground X 2 II state to the A 2 £ + state is calcu- 
lated to be about 4-10 -7 . This extremely small value 
reflects the ionic character of the two states. The 
transition is between different orientations of the 0~ 
ion which has no dipole moment. 
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Table 5. LiO spectral distribution from X 2 I1 and A 2 2+ states 



Transitions 


a 


|<r>| 


|<V>| 


/(r,V) 




Final 
state 


Ion 


Ryd- 

berg 


2 n -* 


( 2 n ; 


3 n, 


2 2+) 


0.0217 


0.00271 


3.94 x 10- 5 


2 n -> 


( 2 n ; 


3 n, 


2 X+) 


.0193 


.00639 


8.26 x 10- 4 


m -> 


( 2 n ; 


'n, 


2 2 + ) 


.125 


.0156 


1.31 x 10- 3 


2 n -* 


( 2 n ; 


i n, 


2 2+) 


.111 


.0368 


2.75 X 10- 3 


2 n -> 


( 2 2 + , 


'n, 


2 n) 


.0415 


.0179 


4.97 X 10- 4 


2 n -> 


( 2 s-, 


'n, 


2 n) 


.0415 


.0179 


4.97 X 10- 4 


m -» 


( 2 2+; 


3 n, 


2 n) 


.0239 


.0103 


1.65 X 10- 4 


2 n -* 


( 2 S-; 


3 n, 


2 n) 


.0239 


.0103 


1.65 X 10- 4 


2 n -> 


( 2 A; 


'n, 


2 n) 


.0587 


.0254 


9.94 X 10- 4 


2 n -> 


( 2 A; 


3 n, 


2 n) 


.0678 


.0293 


1.32 X 10- 3 


2 n -> 


( 2 X~, 


3 S- 


, 2 2+) 


.207 


.0209 


2.89 X 10- 3 


2 n -» 


( 2 S", 


3 S- 


, 2 2+) 


.258 


.0812 


1.40 X 10- 2 


2 n -> 


( 2 A; 


! A, 


2 2+) 


.254 


.0256 


4.34 X 10- 3 


2 n -» 


( 2 A; 


! A, 


2$+) 


.316 


.0995 


2.10 X 10- 3 


2 n -* 


( 2 n ; 


3 £- 


2 n) 


.243 


.0279 


4.52 x 10- 3 


2 n -» 


( 2 n ; 


'A, 


2 n) 


.297 


.0342 


6.79 x 10- 3 


2£+^ 


( 2 n ; 


j n, 


2 S) 


.308 


.0548 


2.26 x lO- 2 


2 S+ ^ 


( 2 n ; 


j n, 


2 S) 


.236 


.1012 


3.18 x 10- 2 


2 S + -^ 


( 2 2 + ; 


3 n, 


22) 


.0535 


.0095 


3.39 x 10- 4 


2 £ + ^ 


( 2 2 + ; 


3 n, 


2J) 


.0409 


.0175 


4.79 x lO 4 


2X+^ 


( 2 £ + ; 


■n, 


2 n) 


.143 


.0111 


1.06 x 10- 3 



Table 5. LiO spectral distribution from X 2 II and A 2 £ + states — 
Continued 



Transitions a 



Final Ion Ryd- 
state berg 



|<r>| 



Kv>l 



/(r, V) 



2 2+-> ( 2 2+; 3 II, 2 n) 0.248 0.0192 3.19 X 10- 3 

2 £+^( 2 II; 'X^II) .0571 .0226 1.72 X 10~ 3 

a The transition 2 n -» ( 2 II; 3 n, 2 2 + ) designates the following: 
initial state — > (final state; ion state, Rydberg orbital symmetry). 

4. References 

[la] Wahl, A. C, (private communication). Energy and spectro- 
scopic constants reported in reference [2]. 
[lb] Yoshimine, M, J. Chem. Phys. to be published. 
[2a] Freund, S. M., Herbst, E., Mariella, R. P., and Klemperer, W., 

J. Chem. Phys. 56, 1467 (1972). 
[2b] White, D., Seshadri, K. S., Dever, D. F., Mann, D. E., and 
Linevsky, M. J., J. Chem. Phys. 39, 2463 (1963). 
[3] Huzinaga, S., J. Chem. Phys. 42, 1293 (1965). 
[4] Liu, B., private communication, to be published. 
[5] Hansen, A. E., Mol. Phys. 13, 425 (1967). 
[6] Schadee, A., J. Quant. Spectr. Rad. Trans. 7, 169 (1967). 



(Paper 76A6-751) 



668 



